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ABSTRACT 

The predictions of the atmospheric v event rates are affected by significant uncer- 
tainties, however the evidence for the 'disappearance' of z/^'s and F^'s obtained 
by SK and other underground detectors is robust and cannot be accounted in 
the framework of the minimum standard model without assuming very large 
ad hoc experimental systematic effects. The existence of 'new physics' beyond 
the standard model is therefore close to be established; v oscillations provide a 
very good fit to all data. The theoretical uncertainties do have an important 
role in the detailed interpretation of the data, and in the estimate of oscillation 
parameters. 



1. Introduction 

The experimental results of Super-Kamiokande and other detectorsi'lSi 

on atmospheric neutrinos have provided strong evidence for the existence of u oscil- 
lations or other forms of 'new phyics' beyond the minimum standard model. This 
conclusion comes from the comparison of the experimental results with theoretical 
predictions that have schematically the structure (for /x-like events as an example): 

,p '^^ = E [<^0 ® ^p.n^^c ® P{^a Z^m)] ® (g) Adet (1) 

where N^^ is the predicted event rate, and VL^ are the muon energy and direction, 
and the different ingredients of the prediction are: (i) 0o the primary cosmic ray (c.r.) 
flux; (ii) Yp(^n)^va the z/ yield, that is the number of neutrinos of flavor a produced 
in p (n) induced c.r. showers; (iii) a^, the neutrino cross sections; (iv) P{Pa ^p) 
the oscillation probability; and (v) Adet the detector acceptance. We have dropped 
all dependences on flavor, energy and zenith angle of the different elements and the 
symbol ® indicates the need for appropriate convolutions that can be done accurately 
only with Montecarlo methods. 

The fundamental result obtained with the measurements of atmospheric u^s is 
that the data cannot be described in the minimum standard model, (in the absence 
of oscillations). In particular the survival probability -P(t'^ — > Vp) must be smaller 
than unity in a large region of energy and v pathlenghth. In short: muon neutrinos 
disappear. 

It is possible to go beyond this fundamental but qualitative statement. With 



the inclusion of //-oscillations one can obtain a good description of the data and 
extract quantitative information on the v masses and mixing" Two flavors v^, ^ 
Ur oscillations with maximal of very large mixing are a viable explanation for the 
data. The determination of the u oscillation parameters is of course of fundamental 
importance, and has also a crucial role to determine the best strategy for future 
experimental studies. 

In this work we will analyse the first three elements of the predictions: the primary 
flux, the modeling of hadronic showers and the u cross sections (without considering 
at all the important questions of detector acceptance and background estimates), and 
discuss what are the effects of the existing uncertainties. We will argue that the fun- 
damental result that the minimum standard model fails to describe the atmospheric 
u data is robust, and that the evidence for the disappearance of z/^'s cannot be "re- 
absorbed" taking into account theoretical uncertainties. On the other hand these 
uncertainties can play a significant role in the interpretation of the data and in the 
determination of the oscillation parameters. 

In this work we will refer to some published calculations of the atmospheric u 
fluxes. Two detailed Monte Carlo (MC) calculations of the fluxes performed by the 
Bartol groupcl and by Honda, Kajita, Kasahara and Midorikawa (HKKM)§ have 
been used to compute predictions for the interpretation of the existing data. Other 
calculations have been performed analytically by Bugaev and NaumovfiJ or again with 
MC methods by Lee and Kohli. A new detailed calculation using the FLUKA code^l 
for the modeling of hadronic interaction is close to completionli^. 

This discussion is organized as follows. In the next section we will briefly discuss 
the evidence for new physics in the atmospheric u data; in section 3 we will describe 
some predictions for the atmospheric u fluxes; in section 4, 5 and 6 we will consider the 
three main elements in the prediction: the primary cosmic ray flux, the development of 
hadronic showers; and the u cross sections; in section 7 we will consider the constraint 
coming from atmospheric muon measurements; in section 8 we discuss discuss the 
question of the absolute normalization of the ly flux; in section 9 we will discuss some 
possible effects of the uncertainties in the predictions for the determination of the 
oscillation parameters. 



2. Evidence for the disappearance of z/^^'s 

The evidence for the disappearance of muon (anti)-neutrinos comes from the 
observation of three experimental effects, listed here in order of 'robustness': (i) the 
detection of an up-down asymmetry for the /i-like events, (ii) the detection of a 

" Other mechanisms, different from oscihations, such as decay, FCNC, or violations of the equiv- 
alence principle, have been proposed as explanations of the disappearance of the muon neutrinos. 
These alternative or 'exoLtic' models provide (at least until now) significantly less satisfactory de- 
scription of the data. SeeQu for more discussion. 



small /i/e ratio, (iii) the detection of a distortion of the zenith angle distribution and 
a suppression of the v induced upward going muon flux. 

Some comments about the first two effects can be made already before entering 
in a detailed discussion of the prediction because two properties of the v fluxes: an 
approximate up-down symmetry, and a strongly constrained v^j v^. ratio, are to a large 
extent independent from the details of the calculation and provide 'self calibration' 
methods. 

(i) The V fluxes are predicted in the no-oscillation hypothesis to be approximately 
up/down symmetric: 

(i?„7r-0) (2) 

This follows as a simple and 'purely geometrical theorem from two assumptions: (a) 
the primary c.r. flux is isotropic, (b) the Earth is spherically symmetric. The geomag- 
netic field, as will be discussed in some detail in section 4.2, is the most important 
source of up-down asymmetry. The theory predicts unambiguosly that the effects 
of the geomagnetic field are small, decreasing with Ey and for Kamioka produce an 
excess of up-going particles. This predictions have received an important quanti- 
tative experimental confirmation with the measurement of the east-west effect in 
SkIiI, that is also produced by the geomagnetic field, and is largely independent from 
oscillations^. For the sub-GeV and multi-GeV muon samples in SK the up/down 
asymmetry A = {U — D)/{U + D) (where U {D) is the number of up (down) going 
events with cos^^ < —0.2 (cos^^ > —0.2) takes the values: 

= -0.150 ±0.028 ±0.01 = -0.311 ±0.043 ±0.01 (3) 

these are deviations of 5.4 and 7.2 standard deviations from no asymmetry. 

(ii) The and Ue fluxes are strictly related to each other because they are produced 
in the chain decay of the same charged mesons (as in 7r+ — »• i^/^Ai^ — >■ t'^ (i^/^J^ee^)): 

(^,^{E,9) = r,{E,,9)x^,^{E,9) (4) 

the factor r^(p){Ej,,d) vary slowly with energy and angle is quite insensitive to the 
details of the calculation. At low energy when essentially all parent muons decay 
before reaching the ground one has rj,(,7) ~ 2. This is in agreement with a naive 
u 'counting' argument. This simple argument is approximately valid because after 
chain decay the three neutrinos produced in the decay of a pion have similar average 
energy (for ultrarelativistic vr"*" — > /^^(z/^z/gC^) decay one has: (E^/Et,) = 0.213, 0.265 
and 0.257 for the z/^ and z/g channels). With increasing energy, because of the 
Lorentz time dilatation, the muons begin to reach the ground before decay and rj,(,7) 
increases, with a stronger (weaker) effect for vertical (horizontal) directions because 
of the differen distance between the primary interaction point and the ground. These 
energy and zenith angle dependence of r (see fig. P are controlled essentially by 
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Fig. 1. pije ijatio of cc-interacting i/'s calculated in the no-oscillation hypothesis using the Bartol □ 
and HKKM N fluxes. 

geometry and reliably calculable. In SK the the double ratio R = {fi/e)data/ {fi'/G)MC 
for the sub-GeV and multi-GeV samples are measured as: 

Rsub = 0.668 to'^oll ± 0-007 ± 0.052 Rmui = 0.663 tomt ± 0-013 ± 0.078 (5) 

where the errors are statistical, systematic and theoretical (as estimated by SK). Com- 
bining in quadrature all errors the significance of the deviations of the double ratios 
from unity are 5.7 and 3.7 standard deviations. The uncertainty in the denominator 
of the double ratios has several sources: the ratio tt^ /tt~ (that determines ratio fe/T^e 
important because of the different cross sections of z/'s and z7's), the contribution of 
kaons^ the shape of the primary spectrum, and the momentum distribution of the 
produced mesons. These uncertainties have estimated0'El@ as of order 5% for a fixed 
value of Ey. Since experimentally the (/i/e) implies an integration over a finite inter- 
val of energy and the convolution of the v specra with detection efficiency curves that 
are fiavor dependent, an important source of uncertainty comes from the knowledge of 
the shape of the energy spectrum, and the energy dependence of the v cross section. 
For this reason the uncertainty on {^Ji/e)Mc is larger than 5% and was estimated by 

''In the chain decay t'/j/^^ ^ v^iv ^Vf.e~) the average energy of the three neutrinos is 

(EvIEk) = 0.477, 0.159 and 0.205, therefore for neutrinos from K decay the naive v 'counting' 
argument is not valid. This is compensated by the fact that charged ad neutral kaons are also a 
source of electron neutrinos via the decays K -Kevg- 



SK as 8% and 12% for the sub-GeV and multi-GeV sample. The estimate of the 
value and uncertainty of the denominator in the double ratios is very important for 
the interpretation of the data. In any case even with a very conservative estimate of 
'theoretical errors' it does not appear possible to explain the detected low values of 
R obtained without appealing to experimental errors in flavor identification. 

(iii) For the interpretation of the upward-going muon data, the methods of comparing 
up and down hemispheres, or the i/^ and z/g fluxes are not available, therefore the 
uncertainties in the predictions are more important. The angular distribution of the 
muons, and the ratio between a high energy and a low energy part of the signal (or 
example stopping and through-going muons), are less sensitive to the model (see 
for a critical discussion), and provide important tools for the interpretation of the 
datalli. Upward going muons are the best way to study the z/^ fluxes at high Ey, 
and this is important to study experimentally the energy dependence of the v flavor 
transitions!. 



3. Neutrino Fluxes calculations 



In this section we will briefly present some results about the properties of the 
calculated neutrino fluxes. For the sake of clarity, we will discus the v event rate that 
(for example for e-like events) can be calculated as: 

Na [MEu. cos^,) <^(i?.) + <hXE,, cose,) ag^iE,)] (6) 



dEu d cos 01, 



as the product of the flux times the neutrino cross sections. For (Tj, we will use in this 
work the estimates of ill including nuclear effects. 

In fig. 1^ we show the energy distributions of up-going and down-going interacting 
neutrinos at Kamioka as calculated for minimum solar activity, by the Bartoli an 
HKKM0 groups. The u direction was integrated over the angular regions: cos 9^ E 
[—1,0] for 'Up' events, or cos^^j^ G [0,1] for 'Down' events. Plotting E^dN/dE,, = 
dN/dlnEu with a logarithmic energy scale, the area below the curve is proportional 
to the total rate. Some remarks about fig. |^: 

(i) The two calculations, especially for E^ > 1 GeV, are quite similar to each other 
both in absolute normalization and in the shape of energy distribution. This simi- 
larity, as it will be discussed later, is to some extent the result of a cancellation of 
a higher (lower) primary cosmic ray flux and a lower (higher) u yield per primary 
particle for the HKKM (Bartol) calculation. 

(ii) Both calculations predict that the up-going and down-going u fluxes become 
equal for E^, > 1 GeV. At low energy (and in the absence of oscillations) the calcu- 
lations predict a higher up-going flux. This is the result of the geomagnetic effects 
and has opposite sign with respect to the effect of oscillations. It should be noticed 
that the Bartol calculation predicts a stronger no-oscillation up-down asymmetry for 



25 




E, (GeV) 

Fig. 2. Energy distributions of up-going (cos^^ < 0) and down-going {cos9i, > 0) cc-interacting 
(ve +T'e)^s and (z^p +X/^i)'s at Kamioka The calculations are for solar minimum and with the Bartola 
and HKKM0 ly fluxes. 

neutrinos with < 1 GeV. 

The calculated u event rates depends on the detector position and the epoch of 
solar activity. This is illustrated in table 1 and 2, that give the charged current 
event rates (with no detector efficiency or requirement for containement included) 
for two detector positions (the Kamioka and Soudan mines) and different levels of 
solar activity. For the Kamioka site, we also compare the two MC calculations of 
Bartol and HKKM. The quantities {Ul, Uh, Dl and Dh) refer to an integration in 
the z/ zenith angle in the regions: cosOi, G [—1,0] for 'Up' events, or cos 6*1, G [0,1] 
for 'Down' events, and an integration in energy < 1 GeV for 'Low' energy and 
Ey>l GeV for 'High' energy. The asymmetry is calculated as A = {U — D)/{U + D). 
Some remarks can be useful: 

(i) For sufficiently large E^ {Ey > 1 GeV) the calculated event rate is: (a) independent 
from energy, (b) independent from the detector position, and (c) up/down symmetric. 
Conversely at low energy there is a significant dependence on detector location and 
solar epoch, and there are significant deviations from up/down symmetry even in the 
absence of oscillations. 

(ii) At the epoch of solar minimum, that to a good approximation applies to the data 
taking period of SK (see fig. 0), without considering the effect of detector efficiency, 
the rates of e-like events with E^ < \ GeV {Ey > 1 GeV) in Soudan and SK are 
predicted by the Bartol calculation to be in the ratio: Soudan/Kamioka ~ 1.72 (1.17). 



Table 1. Event rates (kton yr) ^ and Up/Down asymmetry for e like event. 



Site Flux 


Ul Dl Al 


Ufj Dlr A%r 


Kamioka HKKM min 
Kamioka HKKM mid 
Kamioka HKKM max 


15.8 13.5 0.077 
15.1 13.1 0.070 
14.3 12.7 0.060 


14.5 14.4 0.002 
14.2 14.2 0.001 
13.9 13.9 0.000 


Kamioka Bartol min 
Kamioka Bartol max 


16.2 12.5 0.129 

14.3 12.2 0.077 


14.6 14.2 0.013 
14.3 14.1 0.006 


Soudan Bartol min 
Soudan Bartol max 


20.7 28.4 -0.158 
17.5 21.9 -0.112 


16.3 17.6 -0.039 
15.8 16.8 -0.031 


Table 2. Event rates (kton yr) ^ and Up/Down asymmetry for /i like event. 


Site Flux 






Kamioka HKKM max 
Kamioka HKKM mid 
Kamioka HKKM max 


27.5 24.4 0.059 
26.5 23.8 0.053 
25.2 23.0 0.046 


34.3 34.2 0.001 
33.8 33.8 0.000 
33.1 33.2 0.000 


Kamioka Bartol min 
Kamioka Bartol min 


27.8 22.3 0.109 
25.1 21.9 0.067 


33.9 33.4 0.008 
33.5 33.3 0.004 


Soudan Bartol min 
Soudan Bartol max 


34.1 44.8 -0.136 
29.8 36.5 -0.100 


36.4 38.4 -0.026 
35.8 37.3 -0.021 



The experimental results suggests a lower ratio (see section 7 for more discussion). 

(iii) All event rates should be reduced during the future period of solar maximum, 
with a significantly larger reduction for Soudan because of the lower geomagnetic 
cutoff. As an illustration, the ratio (0 max)/(0 min) for sub-GeV e-like events is 
predicted by the Bartol calculation to be ~ 0.80 for Soudan and ~ 0.93 for Kamioka. 
These time variations should be measurable and offer a handle for a check of the 
predictions. 

(iv) At low energy the geomagnetic effect result in a significant up-down asymmetry 
even in the absence of oscillations. It is remarkable that this asymmetry has opposite 
signs at Soudan (with an excess of down-going z/'s) and Kamioka (excess of up-going 
i/'s). This effect is also illustrated in fig. ||. The difference between Soudan and 
Kamioka can be easily understood qualitatively taking into account the fact that the 
magnetic latitudes of Kamioka and Soudan are Xm — 27° and 56°, and the different 
geomagnetic cutoffs for the two locations (see section 4.2). 

3.1. The response function 

Fig. § is taken from Gaisseiil and shows the energy distribution of the primary 
particles that (would) contribute to the to the /i-like sub-GeV signal in Super- 
Kamiokande. The solid (dotted) curves are for solar minimum (maximum). The 
pair of curves indicated as B (C) are the response curves for up-going (down-going) 




neutrinos at Kamioka. One can notice how the lower flux of the down-going emi- 
sphere originates from the higher geomagnetic cutoffs relevant for the low magnetic 
latitude site of Kamioka. Curve A is calculated assuming zero cutoff and is a reason- 
able approximation for the down-going v flux of the high magnetic latitude of Soudan 
where geomagnetic cutoffs are small. The primary energies that are important for 
the atmospheric v signal are in a broad range Eq ~ 3-100 GeV. It is also important 
to notice that experiments made in different locations will be sensitive to somewhat 
different ranges of primary energy, and that in particular the event rates at Soudan 
receive a large contribution, from particles with E'o — 10 GeV, and the prediction is 
very sensitive to the flux and interaction properties of low energy primaries. 



4. The primary cosmic ray flux 

The primary cosmic ray flux has been a major source of uncertainty (see@ for a de- 
tailed discussion) because of the discrepant results obtained by two groups: Webbeiil 



In more recent time there have been 
CAPRICEli, BESSS and MASS2§) 



and LEAPBa), differing by ~ 50% (see fig. | 
four new measurements of the p flux: IMAX 
of the proton flux (some including also measurement of higher mass primaries). All 
four these new measuremens (shown in fig. ^) are, at least qualitativelvjin agreement 
between each other and confirm the lower normalization of the LEAPna experiment. 
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Fig. 4. From GaisseitJ. Energy distribution of primary particles that produce the sub-GeV fi-like 
event rate in Super-Kamiokande (see text). 



The conclusion is that, if one is ready to discard the early measurement of Webbei", 
the primary flux is much better determined, with an uncertainty in the relevant region 
of better than 10%. 

It is important to compare the new experimental results on the p flux with the 
representations chosen for the calculations of the u fluxes. This comparison is shown 
in fig. ^ and in more detail in fig. |[ where the data and the different descriptions of 
the p flux are shown as a ratio with respect to the proton flux used in the HKKM0 
calculation for minimum solar activity. The three thin solid lines describe the p flux 
assumed for by HKKM for minimum, medium and maximum solar activity (upper 
line with constant value of unity, middle line and lower line). The thick line is the 
description of the p flux at solar minimum used in the Bartol calculation. Only 
two sets of data points are shown, those obtained by the Capriceil and MASS2E3 
experiments. These data were obtained in 1997 and 1991 approximately at the epoch 
of minimum (maximum) solar activity. Some remarks: 

(i) The measurements are consistent with each other for Ep > 10 GeV. Below this 
energy the difference in the measured flux, agrees with the expectations for the solar 
modulation. 

(ii) For Ep > 10 GeV, the descriptions of the p flux used in the calculations of the z/ 
fluxes are in poor agreement with the data. The discrepancy is not a constant factor. 
As an illustration in fig. ^ we have included a dot-dashed line: {Ep/7 GeV)~°'^^, that 
corresponds very roughly to the ratio data/model for the solar minimum p flux of 
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Fig. 5. Measurements of the c.r. proton flux (points) and the representations used as input in the 
Bartol and HKKM calculations (lines). 



HKKM, and a dashed line (£'fc/9.5 GeV)""-^'*^), that corresponds roughly to the ratio 
Bartol/HKKM at solar minimum. 

4.I. Time dependence 

The time variability of the primary cosmic ray flux is caused by the time vary- 
ing structure of the heliosphere, that is on the electro-magnetic fields that fills the 
interplanetary space and are the result of the interaction of the solar wind (an out- 
flow of plasma from the sun) with the interplanetary medium. The topics of the 
structure of the heliosphere, and the phenomenon of the modulation of cosmic rays 
(the two topics are of course strictly related since the c.r. modulations are an im- 
portant tool to study the characteristics and dynamics of the heliosphere) are very 
complex (see il for a review and references) and will not be discussed here. It can 
only be said that because of the presence of the solar wind plasma and the e.m. 
fields associated with it, for cosmic rays the heliosphere is a medium that goes from 
semi-transparent (for rigidity i? ~ 5 GeV) to totally opaque (for kinetic energy of 
few hundred MeV/nucleon). The transparency of the medium depends on the in- 
tensity of the solar wind, that is correlated with the intensity of the solar magnetic 
activity, and varies with an 11 years period (22 years considering also the polarity 
of the magnetic field). This magnetic activity can be monitored for example with 
measurements of the number of sunspots. The correlation between the cosmic ray 




flux and the solar activity is shown in fig. ^ The intensity of the cosmic rays at the 
Earth in the energy region Eq ~ 1-20 GeV can be continuosly studied at ground level 
using neutron monitor detectors. These detectors measure the hadronic component 
of secondary cosmic rays, detecting the neutrons produced in the interactions of the 
secondary particles in the material (lead) that surrounds the detector. The neutrons 
are moderated (with paraffin or polyethylene) and are detected as thermal particles 
using Boron tri-fiuoride proportional counters. The rate measured by a neutron- 
monitor detectors depends on the magnetic latitude of the detector that determines 
the minimum rigidity of the primary particles that can reach the ground level, and 
by the altitude of the detector, that determines the absorption of the secondary par- 
ticles. The top panel of fig. 0. shows the monthly average of the Sunspot Number 
since 1970, the bottompanel shows the variations in the monthly counting rate of 
three neutron monitorgH^. The two thin histograms that shows a stronger variabil- 
ity give the activity measured by two detectors located at Kiel (altitude 54 meters, 
vertical cutoff 2.36 GV) and Moscow (altitude 200 meters, vertical cutoff 2.43 GV). 
The two histograms are nearly superimposed showing excellent correlation, they are 
measuring real variations of the primary c.r. flux. The thick histogram exhibiting 
a smaller variation is the activity measured at Huancayo (Peru) (altitude 3400 m, 
cutoff Rigidity 12.92 GV). It can be seen how the time variability of the cosmic ray 
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Fig. 7. The top panel shows the monthly average of the Sunspot Number. The bottom banel shows 
the relative rate of three neutron-monitor detectors, Kiel, Moscow and Huancayo (see text). 

flux is reduced to less that 5% for this higher rigidity cutoff. 

4-2. Geomagnetic effects 

The geomagnetic fleld has some very important effects on the propagation of cos- 
mic rays, and forbids low rigidity particles from reaching the surface of the Earth. 

The geomagnetic effects depend strongly on the position on the Earth of the observer, 
and on the direction of observation. For a qualitative understanding let us approx- 
imate the Earth magnetic field as a dipole, and consider two observers, the first at 
one magnetic pole, and the second on the magnetic equator. The 'pole observer' 
looking vertically up, can observe particles of all momenta, since in this case the fleld 

— * 

B and the particle momentum p are parallel and there is no bending of the particle 
trajectories; the observer located at the equator when looking for particles arriving 
horizontally from east can detect protons (or particles with charge +e) only if they 
have a momentum larger than approximately 59 GeV. Particles with this charge and 
momentum have an (unstable) circular trajectory that remains at a constant radius 
r — seeing a magnetic fleld of constant intensity — 0.31 Gauss) always orthog- 



onal to their 3-momentum. It is obvious that particles with lower momentum have 
'impossible' trajectories. With a little more work one can calculate for example that 
the 'equatorial observer' looking horizontally west (vertically up) can observe protons 
with momentum j9 > 10 GeV {p ^ 15 GeV). In general for a detector position and 
direction of observation it is possible to compute a 'cutoff rigidity' Rc{x,fl) 

(0 < i?c ~ 59 GV) such that only particles with R> Rc can be measured. 

The Liouville theorem states that the density of points in phase space volume is 
constant. This fundamental theorem has a simple and important application^ to the 
propagation of cosmic rays in a static magnetic field. Since the momentum of the 
particles does not change in absolute value, then the differential flux along a particle 
trajectory is constant, and it follows that if the flux is isotropic at large distances 
from the Earth, then the flux in a small cone around any trajectory is constant and 
independent from the position. More formally if 0oo(p) is the isotropic flux at large 
distances from the Earth, then the flux measured by a detector at a point x in the 
direction Q has the form: 



4>x{P, ^) = 0oo(p) X C(p, fi, x) (7) 

where the quantity ( is either or 1, but never something else. The cosmic ray energy 
spectrum is therefore not deformed by the geomagnetic field and the effect of is is 
to 'remove' a set of momenta from the spectrum. In summary the calculation of the 
geomagnetic effects, using the facts that field B is approximately constant in time 
and the primary flux is isotropic at large distances, can be reduced to the calculation 
of the allowed and forbidden trajectories. 

This problem can be solved analytically for the special case of a volume that is 
entirely filled with an exactly dipolar magnetic field. In this case the function ( has 
the form of a single sharp step01i. Introducing the rigidity R = p/q the solution is: 

aR,Q,x) = Q[R-Rs{r,X,e,if)] (8) 

where B(x) is the Heavyside function {Q{x) = for x < 1, G(x) = 1 for x > 1), r is 
the distance from the dipole center, A is the magnetic latitude, 6 is the zenith angle, 
(fi is the azimuth angle measured clockwise from the magnetic north, and Rs is the 
Stormer0 rigidity cutoff that for positively charged particles is: 

Mr,X,e,^) = {^^) { + j3°A',,t^,in^)i/2]2} (9) 

(for negatively charged particles the cutoff is obtained with the reflection (p ^ ip + Tc) 
that is exchanging east and west). In (^ M is the magnetic dipole moment of field. 
For the Earth M ~ 8.1 x 10^^ Gauss cm that corresponds to a polar magnetic field of 
0.62 Gauss. The quantity M/2r^ ~ 59.4 GV corresponds to the rigidity of a particle 
in a circular orbit of radius in the earth's magnetic equatorial plane. 



The Stormer analytic result is only an approximation for a realistic calculation 
because of two effects: (i) the Earth field is not exactly dipolar, (ii) the field does 
not fill the entire space, and therefore a fraction of the trajectories that are 'allowed' 
according to the Stormer formula intersects the surface of the Earth and is in fact 
forbidden. The problem can be solved in general numerically using the backtracking 
method. In this method to determine if a cosmic ray trajectory {x,p, q) is allowed or 
not, one integrates numerically the equation of motion for a particle of momentum 
—p and charge —q, starting from the point x using a detailed description of the ge- 
omagnetic field, this of course corresponds to the calculation of the past-trajectory 
of the particle considered. If the trajectory 'goes to infinity' without intersecting the 
surface of the earth, the initial trajectory is allowed: CiPi^i^) = 1- If the trajectory 
intersects the Earth, or remains confined within a finite radius the trajectory is for- 
bidden: CiVi^il) = 0. We note that in the exact calculation the function ({R,Q) 
has not anymore the form of a single step function, and in a narrow range of p/q the 
are alternating intervals of allowed and forbidden rigidities (the 'penumbra' region), 
therefore strictly speaking for a fixed position x and direction Q one should not speak 
about a single cutoff rigidity but of a minimum allowed rigidity Ri and a maximum 
forbidden rigidity R2 (with Ri < R2), however since i?i ~ i?9 the concept of cut- 
off remains qualitatively useful. Both the Bartolil and HKKM0 calculations use the 
backtracking method to estimate the geomagnetic effects. An important confirmation 
of the correctness of the estimate of the geomagnetic effects has been obtained with 
the experimental observation of the east-west anisotropy for atmospheric neutrinos 
performed by SK0. 



4-3. Isotropy of the c.r. flux 

In the previous subsection we have assumed that the c.r. flux , in the absence of 
geomagnetic effects is isotropic. This 'intrinsic' isotropy can in principle be studied 
measuring the flux in a fixed direction in the detector system of coordinates so that 
the geomagnetic effects are constant. Since the Earth is rotating this allows to study 
the variations of the primary flux in a cone in celestial coordinates. The measure- 
ments for the range of primary energy of interest have shown the existence of 'solar' 
anisotropies (related to day-night effects) with amplitudes of a fraction of a percent 
and no 'sidereal' anisotropiesEj. 



5. Modeling of the development of hadronic showers 

The modeling of the development of hadronic showers is a second important source 
of uncertainty in the calculation of the u fluxes. A cosmic ray shower has a 'tree' 
structure and can be seen as made of branches and vertices. The primary particle 
interaction point is the first vertex of the tree, the secondary particles produced in 
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particle is destroyed and several 'out' particles (or new branches) are produced. The 
new particles can also be followed, and the process recursively iterated. The neutrinos 
are produced in the weak decays of kaons and charged pions. Because of this tree 
structure, the calculation of cosmic ray shower is a problem that can be solved in the 
most natural and accurate way with Montecarlo techniques. Analytic treatments are 
also possibleinijii and represent a reasonable approximation that allows to gain insight 
on the physics of the problem, however the MC methods are more general and precise 
and have to be preferred. 



5.1. Modeling of the hadronic interactions 

The main uncertainty in the modeling of hadronic showers is the description of 
particle production in hadronic interactions. The multiplicity, energy spectrum and 
flavor composition of the particles produced in the final state are all important factors 
in the determination of the u fluxes. 

In fig. H taken fromS a review paper on the calculation of atmospheric u fluxes 
we show the energy distribution of charged pions produced in the interactions of 
24 GeV protons on air as modeled in different calculation of the atmospheric u fluxes. 
Inspection of fig. ^ shows that the Bartol calculation has a higher multiplicity and a 
slightly harder energy spectrum of the HKKM one, and much larger difference with 
respect to the Bugaev and Naumov estimates. This difference in the modeling of 
particle production would account for a ~ 20% difference in the u event rate between 
the Bartol and HKKM calculation if the same primary flux is assumed. 

A comparison of the FLUKAi and Bartola interaction models, shows that the 
FLUKA code also predicts a smaller multiplicity for charged pions, this will result 
in a smaller u flux for the same primary flux. The details of the full calculation will 
soon be made available. 



5.2. The method of calculation: 3D versus ID 



The MC calculation of an hadronic shower is in principle a straightforward, even 
if technically not trivial, computational problem. All the physics that is needed is 



well known, with the only exception of the description of the hadronic cross sections, 
discussed in the previous subsection. However all the calculations that are publi- 
cally available (with the exception of a first attempt by Lee and Koh0) have been 
made in the approximation of considering the development of the hadronic showers 
as one-dimensional. This implies (discussing in particular the implementation of the 
approximation in the Bartol code) the following: (a) multiple scattering and devia- 
tions in the magnetic field are neglected: all particles propagate along straight lines; 
(b) in the generation of the interaction or decay vertices the kinematical properties 
of final state particle are calculated exactly including the generation of p± and the 
exact conservation of 4-momentum, however after the generation of the vertex is 
completed, for all particles with p\\ > the angle with respect to the primary particle 
is set to zero while particles with py < are disregarded. In this way each u is exactly 
collinear with its primary particle. 

The motivations for this approximation are purely technical: the saving of com- 
puter time. The neutrinos that reach a detector in a certain geographycal position 
(say Super-Kamiokande) are in general produced in cosmic ray showers that devel- 
oped anywhere on the earth and (allowing for z/'s non collinear with the primary 
particle) with any direction. It is clear that the Montecarlo generation in this gen- 
eral scheme is very inefficient, because one has to generate all possible cosmic ray 
showers in the Earth atmosphere, but only a very small fraction of them will produce 
neutrinos in the vicinity of the detector. With the assumption of collinearity only 
a very small subset of the cosmic ray has to be studied and the problem becomes 
numerically enormously simpler. 

The dominant contribution to the angle between the neutrino and the parent 
particle arises because of the non-coUinearity of the parent mesons with the primary 
particle, and (for u from muon decay) because of the deviation of the /i^ is the 
geomagnetic field with the ttu angle a smaller contribution: 

(^^0.) - © ^^(^^ 5e) - ^ © f - © (10) 

A first test of the importance of 3D effects in the calculation of the atmospheric v 
fluxes has been performed with FLUKA0. The test is a comparison of two calcula- 
tions of the z/-fluxes performed with the ID approximation or with a fully 3D method, 
but in both cases assuming spherical symmetry for the Earth and the absence of the 
geomagnetic field. In this approximation, all points on the surface of the Earth are 
equivalent and can be considered as a single detector, therefore the CPU time problem 
is not present, note however that since the deviation of the muons in the geomagnetic 
field is neglected, the angle between v and primary particle is understimated. The 
results of the test show that the effects of the 1-dimensional approximation are neg- 
ligible for multi-GeV events, and a small correction for sub-GeV events. A complete 
3D calculation is feasible using an efficient weighting of the simulated events. Results 



on fully 3D calculations of the v fluxes should very soon be made publically available 
(at least) from the Fluka and Bartol groups. 

5. 3. The density of the atmosphere and the mountains 

The density of the medium where the c.r. shower develops, changes with continuity 
and it is essential to take this fact into account. The density profile of the atmosphere 
can to a good approximation be taken as constant in time and independent of position 
(the concept of altimeter is based on this approximation). Variations of the average 
density profile with the geographical latitude and the season have been investigated 
and are the source only of minor corrections for low energy neutrinos. For accurate 
calculations it is also necessary to consider the profile of the mountains above the 
detector. This is straightforward to do, and represents a few percent reduction of 
the component of high energy and vertical down-going neutrinos that is produced in 
muon decay, the effect is therefore stronger for electron neutrinos. 

6. The neutrino cross section 

The description of the neutrino cross section is a source of uncertainty of compa- 
rable importance to the primary flux and the hadronic cross sections in the prediction 
of the atmospheric v event rates. At high E^, when most of the phase space for u 
interactions is in the deep-inelastic region, (jy is reliably calculable in terms of well 
determined parton distribution functions (PDF's). However in the energy range rel- 
evant for atmospheric neutrinos {E^ ~ 1 GeV) the description of (Ti, is theoretically 
more difficult, and the available data is not sufficient to constrain the cross section 
to better than 15-20%. 

Quasi-elastic scattering is the most important mode, its importance being en- 
hanced by the event selection criterion (the 'single ring' condition) that requires a 
single visible particle in the final state. Uncertainties on this process arise from uncer- 
tainties in the axial form factor of the nucleons (that cannot be measured in electron 
scattering), and more important from the description of the nuclear effect corrections. 
Most MC codes use some version of the Fermi gas model 0. More sophisticated stud- 
ies have attempted to calculate the nuclear effects beyond the Fermi gas moded 
using different theoretical methods to describe the nuclear medium. These studies 
were originally developed to see if nuclear effects could somehow result in different 
cross section for z/g's and z/^'s, and so explain the low (/i/e) ratio observed experimen- 
tally. No such flavor dependence has been found, however the calculations have also 
estimated significant effects (10-15%) in the absolute value of the QE cross sections. 
The Soudan detector in a good fraction of the events can measure together with the 
charged lepton, also the recoil proton in processes as: t'^Fe — > ii^pA*, where A* is 
an undetected nuclear system. The resolutions that can be attained for the neutrino 




kinematical properties {Ei^ and Q^,) depend on how much 4-momentum is carried 
away fron the undetected nuclear system (the 'spectator' in the simple Fermi gas 
model), and can depend on the modeling of the nuclear effects. A detailed study 
should investigate if the resolutions estimated in the Fermi gas modeli are correct. 

Processes where the target nucleon is excited to to a resonance are also of great 
importance in the energy region relevant for atmospheric neutrinos. Both the Super- 
Kamiokande and Soudan MC codes describe this cross section using the algorithms 
of Rein and SehgalS that apply to neutrino interactions the model of baryonic res- 
onances of Feynman, Kislinger and Ravndal0 that describes the baryons as bound 
states of three valence quark in an harmonic oscillator potential. Two examples of 
the prediction for da/dW {W is the mass of the hadronic system in the final state), 
calculated according to the original work, are shown as an illustration in fig. We 
note that the Rein and Sehgal model predicts a complex structure in the W distri- 
bution, and at intermediate energy not only the A resonance (the strongest peak at 
W ~ 1.232 GeV) but also other states give important contribution to the cross sec- 
tion. The existing data does not give a clear confirmation of the prediction. Nuclear 
effects are also important for resonance production. In particular in a nuclear medium 
there is the possibility of a baryonic resonance that is dis-excited to a nucleon, with 
the scattering appearing experimentally as QE scattering. 

All other channels in the u cross section are usually described using the standard 
expression for deep inelastic scattering (DIS)@. In fig. [TO we show and examples of the 
differential cross sections da/dW and daJdlnQ'^ calculated using the DIS formulae, 
and the leading order PDF's of GRV94e3. Looking at da/dW in the left panel of 
fig. ITU] one can see that the as expected the DIS formula does not 'know' anything 
about the hadronic mass spectrum, so it does not describe resonances and in fact it 
computes a cross section also in the unphysical region nip < W < {nip + niT,) that is 
forbidden because no hadronic state with mass in this interval can exist. It remains 
an open problem how to 'join' smoothly the baryon resonance production with the 
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Fig. 10. Left: plot of da/dW for z/^n nX interactions. Right: plot of da/dlnQ^ for the same 
reaction. The ly cross section is calculated using the DIS formulae and the GRV94LO PDF's. 

all other channels described collectively by DIS. It should also be noted, that because 
of the low energy of atmospheric i/'s the in the interactions is always small, even 
if the final state hadronic system is excited above the resonance region (see the right 
panel of fig. |10|) and therefore the use of the DIS framework is not on a solid ground. 

The importance of the neutrino cross section in the prediction of the atmospheric 
neutrino event rates is revealed by the effect of a relatively small modification of the 
description of a,^ in the SK Montecarlo: the choice of a new set of PDF's (GRV94L0i 
replacing the CCFR parametrization) to describe the DISpart of a^. Combined with 
the use of the medium solar activity u fluxes of HKKMliy (that has a small impact 
for multi-GeV events) has led for example to an increase of 4% (8%) for the /i-like 
(e-like) fully contained events, and 7% for the partially contained events (compare 
the MC predictions i lE audi). 

It appears very difficult to calculate accurately from first principles in the 
relevant energy region. The existing data do not determine the absolute value of 
and the energy spectrum of the final state lepton better that ~ 15%. A significant 
improvement can probably only be obtained with the analysis of more data. The 
K2K long-baseline u beanB with a spectrum not too different from the atmospheric 
one offers very interesting possibilities. The placement in the 'near' position along 
the beam of detectors with higher resolution and/or an iron or argon target could 
also be useful for a more complete study. 

One final remark on the question of the u cross section, is that all the experimental 
groups that have made measurements of atmospheric z/'s have compared their results 
to predictions calculated using different descriptions of and different MC codes (in 
general not publically available) for the generation of the u interactions. A detailed 
comparison of the descriptions of ai, used by the different groups is missing. This 
lack of information can be a source of confusion. For example it is important to 
understand the reasons for the difference in the ratios Cdata/eMC stated by the SK 



and Soudan (see section 8) and establish it is due to the v flux or, at least in part, to 
different assumptions about the z/cross sections. 



7. The muon flux constraints 

The fluxes of atmospheric muons are strictly related to the neutrino ones, because 
almost all z/'s are produced either in association with, or in the decay of /i^. Precise 
measurements of the muons offer then a way to check directly the v fluxes. To 
constraint the flux of neutrinos with energy Ey^ the most interesting muons are those 
created with energy around 2 E^, (the factor of two being a kinematical effect). 
The most important range of u energy is E^ ~ 0.1-10 GeV, and especially for the 
lowest part of this range, only a vanishing or small fraction of the related muons reach 
the ground because of energy loss (the atmosphere corresponds to a loss of ~ 2 GeV 
for a vertical muon) and decay (the decay length is £^ ~ 6.2 p^(GeV) km). For low 
energy muons the measurement is possible only at high altitude. 

The authors of the existing atmospheric u calculations checked that their predic- 
tions for the muon flux (that is 'automatically' calculated together with the u flux) 
with the existing muon dataMS at ground level. To a large extent the agreement 
between the HKKM and Bartol calculations of the u fluxes, that is the effect of a 
cancellation between different estimates of the primary flux and different modelings 
of the hadronic interaction, is not casual, but the result of the /i constraint. Recently 
new measurement silHli^l of the ground level muon flux, have been performed with 
detectors developed for balloon flights. The new data suggests that previous esti- 



mates could be too high, (see for example fig. |Tl]). This is an important point that 
need to be studied in detail. New precision measurements of muons at ground level 
in sites with different geomagnetic cutoffs and at different altitudes would certainly 
be very valuable. 

In recent years, the balloon detectors have performed measurements of the muon 
spectrum during their ascent in the atmospher^i'SO'El. The analysis and discussion 
of these results is still in progress. Preliminary comparisons of these results with MC 
calculations have shown some significant discrepancies (especially for E^ < 1 GeV) 
that require a more detailed study. These are difficult measurements, with the altitude 
of the detector continuosly varying, and the statistical errors are still large because 
of the short time available during balloon ascent, some discrepancies between the 
results of different groups do exist, however these and future measurements at high 
altitude (with detectors on balloons or possibly airplaneB) have a good potential for 
reducing the uncertainties in the flux of atmospheric neutrinos. 

8. The 'Normalization problem' 



In the analysis of their data the SK and Soudan collaboration consider the absolute 




normalization of the calculated v fluxes as a free parameter. In the SK analysis of 
the data in term of ^ oscillations, the v fluxes are fltted to the form 

cP,^{E,,e,) = (l+a) ^'^f{E,,9,) {P,^^,^{E,,9,- sm'29,Am^)) (11) 
^.^{E,,9^) = (l+a)0*f(E.,^,) 

where (p^^^^^^-j are the published no-oscillation fluxes of HKKM (or Bartol) and a is an 
unconstrained normalization parameter. In the analysis of 535 days live days of SK, 
for the best flt point, the normalization parameter had the value a = +0.158, also the 
ratios Cdata/eMC for the sub-GeV and multi-GeV samples, had the values 1.18 ±0.03 
and 1.23 ± 0.07 (statistical errors only). In other words according to the analysis inB 
the the calculated u fluxes of Honda et al.0 are too low by ~ 15-20%. This result 
is within the stated uncertainty in the absolute normalization of the calculation, 
however since the primary c.r. flux used in the HKKM calculation is likely to be 
too large, this relatively high measured neutrino rate has generated discussions and 
speculations, at least in part motivated by the observation, that under the hypothesis 
of Ufj^ ^ Ur oscillations the flux if electron neutrinos is not affected by oscillations, 
however is z/^ ^ transitions exist, it can be enhanced. 

In the analysis of 735 days presented at this conference the SK collaboration 
estimates a lower normalization parameter a = 0.084. This can also be seen observing 
that the ratios Cdata/eMC take the new values 1.06 ±0.03 and 1.08 ±0.05 (for sub-GeV 
and multi-GeV). The signiflcant differences with respect to the previous analysis are 



due to a combination of changes in the MC predictions and of statistical fluctuations 
in the data? The MC prediction has increased because of two effects: (i) the data 
is compared to the v fluxes of HKKM calculated for minimum solar activity, while 
in the 545 days analysiJil the data was compared to medium solar activityf (ii) the 
description of the v cross section has been modified as discussed in section 6 with a 
resulting enhancement of the event rate. The combined effects of these changes are 
an enhancement of the MC rate of 4% (8%) for the sub (multi)-GeV sample. 

The result of the Soudan-2 experiment! for the 'shower' events is Cdata/eMC = 
0.81 ±0.09, for all contained events, (0.86±0.11 for the 'high resolution sample') that 
is an effect of opposite sign with respect to SK. In this case the MC calculation is 
based on the Bartoli v fluxes. The difference in Cdata/ cmc between SK and Soudan is 
a hint of a possible problem, and should be carefully investigated. The neutrino flux 
at Soudan is softer that at Kamioka (see fig. H) and is produced on average by lower 
energy primaries (see fig. ^), therefore the difference between the estimated ratios 
Cdata/eMC cit Kamioka and Soudan is not necessarily an experimental discrepancy, 
but it could be for example an indication that the calculated v flux is too soft. In 
section 8, we will discuss the possible effects of a modification of the shape of the v 
energy spectra. 

9. Determination of the oscillation parameters 

In this section I will try to illustrate (only qualitatively) possible effects of the 
uncertainties in the predictions of the v event rates in the the determination of the 
oscillation parameters in the hypothesis of the existence of ^ oscillations. 

9.1. The Soudan experiment 

The Soudan experimei has measured for the 'high resolution sample' a double 
ratio R = 0.59 ± 0.09. This value is obtained as: 0.59 = 0.74/1.25. It is important 
to note that the denominator of the double ratio (/u/e)MC is significantly smaller 
that the naive expectation of approximately twof The low value of (/i/e)Mc can be 
easily understood as the consequence of different detector efficiencies for //-like and 

The difference in the estinmtcs of the rate of sub-GeV e-hke events for the 535 and 736 days analysis 
has generated speculations^ that the i/ fluxes have a time varying component. The analysis ina 
does not take into account small modifications of the pattern reconstruction program that affect the 
entire integrated data sample. Taking this into account, the evidence for a transient component in 
the v flux becomes much weaker. 

The use of the solar minimum fluxes is certainly well motivated considering the solar activity during 
the data taking period (see fig. |^), the higher p flux implied by this assumption could however be in 
worst agreement with the most recent data (see fig. H) on the primary c.r. flux. 
■^The same argument is also valid for the single track/shower events, where the double ratio is 
R = 0.66 ± 0.11, with 0.66 = 0.70/0.95 



e-like events. Muons are more penetrating than electrons, and the containement 
requirement results in a reduced efficiency for higher energy muons, and therefore in 
a low value for for the /i/e ratio. This implies that the (/i/e) ratio is sensitive to 
the shape of the v energy spectrum with a softer (harder) spectrum resulting in a 
higher (lower) ratio. If the calculated v flux (as an example) falls more steeply with 
energy than the true spectrum, the denominator of the double ratio is overestimated, 
therefore R is underestimated and the estimate of the oscillation parameters is biased: 
for a fixed value of sin^ 26', the range of allowed |Am^| is overestimated. 

As discussed in section 7, there is some experimental indication that the shape 
of the z/ energy spectrum is not predicted correctly because the measured event rate 
for sub-GeV e-like events in Soudan is 15-20% lower than the prediction, while in 
Super-Kamiokande, where {Ei^) is higher because of different geomagnetic effects, the 
detected rate is 6-8% higher than the prediction. If this is interpreted as evidence that 
the V spectrum is flatter than the prediction, the allowed region for |Am^| estimated 
by Soudan should be revised to lower values. 

The Soudan experiment, differently from SK, has not observed a clear zenith 
angle modulation of the /x-like event rate. This can be interpreted as an indication 
that the v oscillation length is sufficiently short, so that also down-going are 
significantly suppressed by oscillations, and can be used to exclude low values of 
Am^. This conclusion is also somewhat model dependent, because for the low energy 
i/'s detected by Soudan, the geomagnetic effects induce a non negligible up-down 
asymmetry, with an excess of up-down events. The 'geomagnetic asymmetry' has to 
be subtracted from the data in order to put in evidence the possible suppression of 
up-going neutrino induced by oscillations. From table 2 we can see for example that 
the Bartol prediction for the 'geomagnetic asymmetry' at Soudan depends on the 
solar cycle epoch and decreases for the flatter spectrum of solar maximum activity. 
Therefore a flatter neutrino energy spectrumll would predict of a smaller asymmetry 
induced by geomagnetic effect, then the experimentally detected asymmetry (that 
is of order A = —0.25 ± 0.16) would be interpreted as more significant, and again 
the allowed interval for |Am^| for the ^ hypothesis would be shifted to lower 
values. 

9.2. The Super-Kamiokande experiment 

In the case of the Super-Kamiokande experiment, the up-down asymmetry rep- 
resents really a 'smoking gun' evidence for new physics. The theory predicts unam- 
biguously that at Kamioka the up-down asymmetry is positive (excess of up-going 
particles) and of order A ~ 0.01 for the multi-GeV events. The detected asymmetry 
becomes more pronounced with increasing energy, precisely as expected considering 
the improving correlation between neutrino and muon direction, and the shape of the 
distortions of the zenith angle distributions also match very well the expectations. As 




a consequence, the possibilities of a very large difference in acceptance for up-going 
and down-going particles, or of the presence of a very large non-neutrino background 
appear extremely contrived and unlikely. 

The up-down asymmetry of multi-GeV events is an excellent method to mea- 
sure the mixing parameter sin^ 26. Since the pathlength for neutrinos coming from 
vertically up or down differ by three orders of magnitude, it can be expected that 
for a broad range of |Am^|, if the correlation between the and the detected 
directions is good, one will have that oscillations are absent for down-going particles 
and can be averaged for up-going particles. Selecting two cones around the vertical, 
it is natural to expect that t/ ~ f/o (1 — sin^ 29/2), and D ~ Z^q — Uq. The mixing 
parameter can then be determined independently from the squared mass difference 
as sin^2^ = 2(1 — U/D) with a precision that is essentially limited by the statistical 
errors. Using this 'algorithm' the SK result for the asymmetry for the multi-GeV 
data: A = —0.31 ± 0.04 can be translated in a measurement of the mixing param- 
eter sin^ 29 = —AA/{\ ~ A) = 0.95 ± 0.09, a result in good agreement with what is 
obtained in the detailed fit. 

This argument however tells us that the up-down asymmetry is not a good quan- 
tity for the measurement of Am?. This is illustrated in fig. 12 that shows as a function 
of |Am^| the number of standard deviations of the asymmetry measurement for the 
sub-GeV and multi-GeV sample from predictions calculated assuming the presence 
of ^ Vj. oscillations with maximal mixing. The shapes of the curves show that 
there is a broad interval of |Am^| where the asymmetries have the measured value. 




In fig. we show as a function of |Am^| the deviations of the (/i/e) double ratios 
for the SK sub-GeV and multi-GeV data with respect to the prediction made assum- 
ing again the presence of z/^j ^ Vr oscillations withn maximal mixing. To compute 
the number of standard deviations the value of ctr has been estimated combining in 
quadrature the statistical errors with the systematic uncertainty as given by SK. The 
shape of the curves in fig. are qualitatively easy to understand, the suppression of 
the /i-like events grows monotonically from zero for very small jAm^l when no neu- 
trinos has time to oscillate, to 0.5 at large |Am^|, when the effect of rapid oscillations 
can be averaged. Since the oscillation probability is a function of |Am^|/i?,^, and 
the measured double ratios for the sub-GeV and multi-GeV are approximately equal 
(0.67 and 0.66), the curve for the multi-GeV data that describes the suppression of 



higher energy neutrinos, is shifted to higher values of of |Am^|. Comparing fig. 12 
and fig |l^ one can see that the lowest value of | Am^| of an allowed interval for SK is 
determined by the measurement of the double ratio for multi-GeV events. 

As in the case of Soudan, the analysis of the double ratio is model dependent, 
because since the e-like and /i-like events are produced by different ranges of E^, dif- 
ferent assumptions for the shape of the neutrino energy spectrum result in different 
predictions for the (/x/e) ratio, and therefore to different interpretations of the experi- 
mental results. In fig. |l^ shows the approximate energy distributions of the neutrinos 
that produce different event samples in SK. It is instructive to observe that because 
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Fig. 14. Approximate energy distributions of interacting neutrinos for different classes of events in 
SK. Note the difference in the shape of the response for /i-hke and e~hke events. 



of the containement requirement, in the the multi-GeV sample fully contained (FC) 
e-like events have a broader energy distribution (extending to higher E,^) than FC 
yU-like events, and in fact the MC double ratio for these events is 1.16 (significantly 
smaller than two). However /z-like events are also detected as partially contained 
(PC) events. For this class of events the single ring requirement that is an important 
source of inefficiency for large E^, when the average multiplicity grows, is removed. 
For this reason the PC events are produced by a still broader range of E^ that extends 
beyond the energy interval for e-like events. 

This question is illustrated with a 'pedagogical example' in fig. |I5[ In this figure 
we have performed three Montecarlo simulations of the SK data that differ only 
for the description of the neutrino spectrum. One calculation uses the standard solar 
minimum HKKM fluxes, while the other two assume a flattening (steepening) of both 
z/e and i>fj_ fluxes with a large distortion factor oc E"^^'^^ (the absolute normalization 
has no importance in this analysis). The spectral distortion has a neglible effect on 
the (/i/e) ratio at a fixed E^, but has a signiflcant impact on the predicted ratio for 
the rates for the entire multi-GeV sample because of the different energy dependence 
of the detection efficiency for the two fiavors. The fiatter spectrum predicts a larger 
no-oscillation (/i/e) ratio, and in the interpretation of the data in terms of z/^ ^ Uj. 
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Fig. 15. Deviation (in standard deviations) of the measured SK multi-GeV double-ratio with 
respect to an expectation calculated assuming the presence of maximal mixing <-> Ur- The 
dashed and dot-dashed curve are calculated assuming a distortion of the HKKM u fluxes. 



oscillations results in a larger |Am | to explain the larger suppression:? 
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10. Summary 

The evidence for the disappearance of i/^'s in in the atmospheric neutrino data 
is robust, and represents a very important result. The detailed interpretation of the 
experimental results is limited by systematic uncertainties in the prediction of the 

event rates. 

One of the main sources of uncertainty related to the spectrum of primary cosmic 
rays, has been significantly reduced thanks to new measurements with detectors on 
balloons. One more measurement of the primary fluxes of the Anti Matter Spec- 
trometer (AMS) should soon become available. This new information on the primary 
fluxes has not yet been included in the simulations. 

Uncertainties in the modeling of particle production in hadronic interactions is 
now probably the most important limitation in the calculation of the atmospheric 
fluxes. New measurements of the properties of particle production in proton- 
nucleus interactions in the energy range Eiab ~ 3-100 GeV could signiflcantly reduce 
the systematic error in the prediction of the normalization and energy spectrum of 



•''An additional efii'ect is also that for a flatter spectrum, the average energy of the neutrinos that pro- 
duce the events in the sample is higher and therefore even to produce the same average suppression, 
the |Am^| must be larger. 



the atmospheric u-Huxes. The improvement in the determination of the primary 
flux makes this experimental program more attractive and useful because a better 
determination of the hadronic cross sections can now be directly translated in an 
improvement in the determination of the atmospheric z/-fluxes. 

The neutrino cross section has also to be recognized as an important source of 
systematic uncertainty, also relevant for the LBL disappearance programs. The best 
perspectives for improvement come again from the performance of new experiments. 

Measurements of the muon flux as a constraint for the u flux calculation are also 
very important. The results of several measurements at high altitude with detectors 
on balloons are in the process of being analysed. New precise measurements of the 
muon flux at ground level can also be very important. 
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